A stochastic rainbow caustic observed with cold atoms 
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We report the direct observation of a novel type of rain- 
bow caustic. In contrast to known examples, this caustic 
originates from a dissipative, stochastic process. We have ob- 
served this using cold 87 Rb atoms bouncing inelastically on 
an evanescent- wave atom mirror. The caustic appears as a 
sharp peak at the lower edge of the asymmetric velocity dis- 
tribution of the bouncing atoms. The stochastic process is a 
spontaneous Raman transition due to photon scattering dur- 
ing the bounce. The results are in good agreement with a 
classical calculation. 
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Caustics are ubiquitous phenomena in nature. Exam- 
ples are the cusp-shaped patterns of light reflection on the 
inside of a coffee-cup and the patterns of bright lines ob- 
served on the bottom of a swimming pool @j. The prime 
example of a caustic is the common rainbow, which can 
be understood in a ray-optics picture by considering how 
the scattering angle of a light ray depends on its impact 
parameter on a water droplet Q . Whereas the incident 
rays have smoothly distributed impact parameters, the 
outgoing rays pile up where the scattering angle has a 
local extremum. Such a divergence of the ray density, 
the caustic, appears at the rainbow angle. In atomic || 
and nuclear ^ ] scattering experiments analogous rainbow 
phenomena have also been observed. 

The examples of caustics that have been known so far 
have in common that the outgoing parameter (scattering 
angle) is a deterministic function of the incoming param- 
eter (impact parameter) . In this Letter we report on our 
observation of a new type of rainbow caustic existing by 
virtue of a stochastic process, which distributes a single- 
valued "impact parameter" over a range of "scattering 
angles" . To our knowledge, such stochastic caustics have 
not been observed before. 

We have observed this caustic in the vertical velocity 
distribution of cold atoms, after bouncing inelastically off 
an evanescent- wave mirror Previous experiments 

on the transverse velocity distribution of atoms bouncing 
elastically on corrugated mirrors ||,[l0| also allow an in- 
terpretation in terms of caustics. However, those were of 
the usual deterministic kind, where the outgoing trans- 
verse velocity is a deterministic function of the position 
where the atom hits the mirror. The caustic then origi- 
nates from atoms reflecting from inflection points on the 



mirror surface. 

Whereas other experiments on inelastically bouncing 
atoms concentrated on the cooling properties, the ap- 
pearence of the velocity caustic seems to have escaped 
attention. The incident atoms in our experiment are 
nearly monochromatic, i.e. they have a narrow velocity 
distribution Av/v = 0.03. During the bounce the atoms 
are optically pumped to a different hyperfine ground- 
state, by a spontaneous Raman transition. They leave 
the surface with less kinetic energy, due to the difference 
in optical potential. The outgoing atoms have a broad 
velocity distribution. This is due to the stochastic nature 
of the spontaneous Raman process, so that atoms make 
the transition at different depths in the evanescent wave. 
Our measurements and analysis show that the resulting 
velocity distribution is highly asymmetric. The caustic 
appears as a sharp peak at the minimum possible veloc- 
ity. This divergence demonstrates the strong preference 
for making the Raman transition at the turning point. 

The experiment is performed in a rubidium vapor cell. 
We trap about 10 7 atoms of 87 Rb in a magneto-optical 
trap (MOT) and subsequently cool them in optical mo- 
lasses to 15 /iK corresponding to a r.m.s. velocity spread 
of <r v = 3.8 cm/s. The MOT has a 1/e 2 radius of 0.6 mm 
and is centered 5.8 mm above the horizontal surface of a 
right-angle BK7 prism. When the falling cloud arrives at 
the surface, the mean velocity is 34 cm/s, and the veloc- 
ity spread is a v = 0.9 cm/s. Note that this distribution 
is for fixed z = 0, not for fixed time. At the surface, we 
create an evanescent wave (EW) by a Gaussian shaped 
laser beam of 0.8 mm 1/e 2 radius, which undergoes total 
internal reflection. 

The blue detuned EW induces a repulsive optical 
dipole potential Up{z) — Up(0) exp(— 2nz), where the 
subscript F = 1,2 denotes the hyperfine ground state, 
k = ko \/ n 2 sin 2 9 — 1 , with ko = 2ir/\ the free space 
wave vector of the light, n = 1.51 the refractive in- 
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mum potential Uf(0) cx Iq/Sf, where Iq is the incident 
intensity at the glass surface, and Sf is the detuning. 
Here we define the detuning Si 2 relative to the transi- 
tion 5S 1/2 (F = 1,2) -> 5P 3/2 (F' = F + 1) of the D 2 
line (780 ran), see Fig.|Tj. Note that we neglect the van 
der Waals surface attraction. This is justified as long 
as Uf{0) is much larger than the kinetic energy of the 
atoms. The main influence of this force is a reduction of 
the effective mirror area Hi- 
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An atom entering the EW in the F — 1 state is slowed 
down by the potential Ui(z). Spontaneous Raman- 
scattering can transfer the atom to the higher hyperfine 
ground state (F = 2). When transferred into this state, 
the potential acting on the atom is 1/2(2), which is weaker 
due to the increase of the detuning by approximately the 
hyperfine splitting 82 ~ Si +Shf s (FigJ^). The ratio of the 
two potentials, (3 = 1/2(0) /£/i(0), quantifies the reduction 
in potential energy. As a result the atoms will bounce in- 
elastically, i.e. to a lower height than their initial release 
position. 

Experimental data of bouncing atoms are taken by ab- 
sorption imaging. After a variable time delay the atomic 
cloud is exposed to 20 /xs of probe light, resonant with 
the 5Si /2 (-F = 2) -> 5P 3/2 (F' = 3) transition. Thus only 
atoms which have been transferred to F = 2 contribute 
to the signal. The atomic cloud is imaged on a digital 
frame-transfer CCD camera. The presented images show 
an area of 2.2 x 4.5 mm 2 with a pixel resolution of 15 /jm 
(Fig.||) . The initial position of the MOT is outside the 
field of view. The horizontal line at the bottom shows 
the prism surface. 

A typical series with 2 ms time increments is shown 
in Fig.|2|, where t = is defined as the time that the 
center of the cloud reaches the mirror. Each density pro- 
file has been converted into a horizontal line sum. The 
solid curve is the result of a calculation described be- 
low. The amplitude of the experimental curves is rescaled 
such that the maximum optical density of the experimen- 
tal curve coincides with the theoretical maximum value. 
This is the only fit parameter. 

As expected, the atoms bounce up less high than 
their initial height. Furthermore, the spatial distribution 
shows another striking feature: it displays a high den- 
sity peak at low z, and a long low-density tail extending 
upward. Note also that there is a time-focus: the den- 
sity peak is sharpest when the slowest atoms reach their 
upper turning point. This density peak is an immedi- 
ate consequence of a caustic appearing at the minimum 
possible velocity. 

The caustic can be understood by considering the 
atoms as point particles moving in the EW potential. 
This corresponds to the ray-optics limit for the optical 
rainbow. We consider an atom arriving at the surface in 
its F — 1 hyperfine ground state, with an initial down- 
ward velocity < 0. Its trajectory through phase space 
is determined by energy conservation: f/i(0) exp(— 2kz) + 
±77w 2 = ^mvf, and is depicted in Fig.||by the thick line. 
While it is slowed down by the EW it may scatter a pho- 
ton at a velocity v p , and be transferred to the F = 2 
state. The atom continues on a new trajectory deter- 
mined by 1/2(2). To illustrate the formation of the caus- 
tic, possible trajectories starting at various positions in 
phase space are depicted as thin gray (v p < 0) and black 
(v p > 0) curves. For asymptotically large z the density 
of curves represents the outgoing velocity distribution, 
showing the caustic where the trajectories pile up. This 
distribution is similar to a rainbow. The density of outgo- 



ing trajectories diverges at the "rainbow velocity". Below 
this velocity the intensity is zero, similar to Alexander's 
dark band in the optical rainbow ||] . Above the rainbow 
velocity the intensity distribution decreases smoothly. 

Despite the similarities, there is a crucial difference 
between a rainbow created by sunlight refracted by wa- 
tcrdroplcts and our 'velocity caustic'. The appearance of 
a rainbow is due to a deterministic process, where the 
scattering angle is uniquely determined by the impact 
parameter. In our experiment the incoming velocity Vi 
is nearly single-valued, and is distributed over a broad 
range of outgoing velocities by the stochastic process of 
spontaneous Raman scattering. 

Given the novel character of this caustic, the analogy 
with known examples can only be valid to some extent. 
For example, the position of the caustic is independent of 
the size parameter k. Similarly the optical rainbow angle 
does not depend on the droplet size. We measured for 
three different values of the decay length kT 1 the atomic 
density profile at the upper turning point (Fig.^Ja). This 
is defined as the highest position reached by the peak 
density. The shape of the cloud did not change, but the 
number of atoms decreased with decreasing k _1 since the 
total number of scattered photons is lower. 

The limitation to the analogy becomes apparent 
through the parameter [3 which determines the position 
of the caustic. This parameter characterizes the "degree 
of dissipation" and therefore has no analogy in the op- 
tical rainbow, or any other deterministic caustic. For 
example, the optical rainbow angle is determined by the 
index of refraction n of water. However, the energy of 
the photons is not changed, and therefore n cannot be 
compared to (3. The change in the position of the caustic 
with P can clearly be seen in Fig]j|b,c. Here we measured 
the height of the upper turning point, as a fraction of the 
initial MOT height, for three values of the detuning. 

In order to quantitatively analyze our experimental 
data, we write the optical hyperfinc-pumping rate dur- 
ing the bounce as T'(z) = (1 — q)TUi(z) /Ti8\, where q 
is the branching ratio to F = 1. We define r](v) as the 
survival probability for the atom to reach the velocity 
v without undergoing optical pumping. This function 
decreases monotonically as 77 = —T'rj, with rj(vi) = 1. 
For \v\ < the solution is r](v) — exp(— (v — Vi)/v c ) , 
where v c = 2kMi/(1 — q)mT, with m the atomic mass. 
When the pumping process takes place at a certain veloc- 
ity v = v p , the atom leaves the surface with velocity Vf = 

\Jvp(l — f3) + ~/3vf (Fig.|l|). This results in a distribution 

of bouncing velocities, resulting from atoms which were 
pumped while moving toward (v~) or away from (v^) the 

surface: n(v f ) = v(Vp (v f ))x\^j\+f](v+(v f ))x\^-\ for 

v^9|wi| <Vf< \vi\ , which diverges at Vf = y/J}\vi\. This 
velocity caustic originates from atoms which are pumped 
near the turning point. There are two reasons for scat- 
tering preferentially at this position. An atom spends a 
relatively long time at the turning point since its velocity 
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is lowest there. In addition, the intensity of the EW, and 
thus the photon scattering rate is highest at the turn- 
ing point. The divergence is an artefact of the ray-optics 
description. It disappears due to diffraction when the 
atoms are treated as matter waves. 

To compare the experimentally obtained spatial dis- 
tributions with the model, we first calculated the one- 
dimensional phase-space density p(z, v). The spatial 
distribution is obtained by projecting p on the z-axis. 
Initially the MOT is described by a normalized gaus- 
sian, p(z, v) cx exp(— {z — z Q ) 2 /2a 2 ) exp(— v 2 /2a 2 ) with 
zq = 5.8 mm, o~ z — 0.3 mm and o~ v — 3.8 cm/s. The 
cloud falls due to gravity and expands due to thermal 
motion. When it arrives at the EW, its velocity distri- 
bution is nearly gaussian centered around v = 34 cm/s, 
with a spread of a v = 0.9 cm/s. When the cloud is re- 
flected, its velocity distribution is changed in the way 
described above. We measure the spatial distribution of 
inelastically bounced atoms after a time of flight t. To 
illustrate the agreement with the experimental data, the 
result for a time t = 7 ms is drawn as the thick line in 
Fig. | 

A very intriguing aspect of our experiment is the pos- 
sibility to observe supernumeraries: interference of two 
trajectories with the same outgoing velocity, but with 
opposite pumping velocity v p . In the velocity distribu- 
tion after the inelastic bounce, the high velocity tail is 
the sum of two contributions: from atoms moving to- 
wards and moving away from the surface when they are 
pumped to F — 2. Interference between these two paths 
should result in oscillations in the velocity distribution. 
This is a nontrivial effect because it involves the sponta- 
neous emission of a photon. 

Given our experimental parameters, we expect a typi- 
cal oscillation period of 1 cm/s. This spacing depends on 
k, just as in a rainbow the supernumeraries depend on 
the droplet size. These oscillations are not yet resolved in 
our measurements, because in the EW-field the magnetic 
substates are not degenerate, each producing a different 
velocity distribution. 

In conclusion, we have observed a new type of caus- 
tic due to the stochastic distribution of a monochro- 
matic input. The caustic appears as a sharp peak at the 
lower edge of the velocity distribution of rubidium atoms, 
bouncing inelastically on an evanescent-wave atom mir- 
ror. The caustic bears some resemblance to the com- 
mon optical rainbow. However its dissipative character 
makes it uncomparable with any deterministic caustics 
and qualifies it as a new physical phenomenon. 
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FIG. 1. (a) Relevant atomic levels of Rb (b) An atom 
enters the evanescent wave in its F = 1 state with initial 
velocity Vi. It is decelerated and spontaneously scatters a 
photon at a velocity v p . After being pumped to F = 2 it 
accelerates and leaves the potential with asymptotic velocity 

Vf. 
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FIG. 2. Absorption images (2.2 x 4.5mm 2 ) at different 
moments after bouncing on an EW, with detuning 92V, and 
ft -1 = 1/im, where t = is denned as the time that the 
center of the cloud reaches the mirror. A line sum shows the 
atomic density distribution above the prism. The dotted line 
indicates the prism surface. The initial MOT was located at 
5.8 mm above the surface. The thick curve in the middle is 
the result of our calculation at 7 ms. 
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FIG. 4. a) Density profiles at t — 14 ms after the bounce 
for k' 1 = 1.96/zm, 0.80^m and 0.62^m, Sew = 150r. The 
curves have been offset by 0.015. Variation of k in this 
range, has no influence on the position and shape of the 
spatial density profile, b) Density profiles for various detun- 
ings: &ew = 67, 150 and 233F, again with 0.015 offset. The 
curves are measured at the moment that the density peak 
reaches its highest point and have the same scale, c) Mea- 
sured height of the upper turning point as a fraction of the 
initial (MOT) height. This represents /3, the relative optical 
potential strength of the F = 1 and F = 2 ground states. 
The curve results from a calculation including the full level 
structure. 



FIG. 3. Construction of the velocity caustic in terms of 
phase space trajectories. The velocity of the incident atoms 
is Vi = —0.34 m/s. The thick solid line shows the trajec- 
tory of the lower hyperfine state bouncing elastically. If the 
atom is pumped to the other hyperfine state it continues on 
a different trajectory. The thin lines represent possible out- 
going paths, each starting at a different pumping velocity, 
depending on the position of Raman-transfer. The density of 
outgoing trajectories diverges, yielding a caustic in the veloc- 
ity distribution. Shown in the upper curve are the total dis- 
tribution (solid line), and the contribution of atoms moving 
towards (dotted) and away from (dashed) the surface while 
being transferred. 
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